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Summary. This paper summarizes results and some open problems about the large- 
scale and long-time behavior of asymmetric, disordered exclusion and zero-range 
processes. These processes have randomly chosen jump rates at the sites of the un- 
(^ . der lying lattice Z'^. The interesting feature is that for suitably distributed random 

^ i rates there is a phase transition where the process behaves differently at high and 

1-^ I low densities. Some of this distinction is visible on the hydrodynamic scale. But to 

c^ ' fully understand the phase transition, results on a finer scale are needed. 



X 



Mathematics Subject Classification: Primary 60K35, Secondary 82C22 
Keywords: Exclusion process, zero-range process, hydrodynamic limit, quenched 
disorder, phase transition 
Short Title: Hydrodynamics for disordered particle systems 



Research partially supported by NSF grant DMS-9801085. 

1 



1. Introduction 

This paper introduces some recent results on the hydrodynamics of disordered, 
asymmetric simple exclusion processes (SEP) and zero-range processes (ZRP). The 
disorder refers to the rates of jumping attached to the sites of the underlying lat- 
tice: The particles move on Z'^, and each site a; G Z*^ has a random variable a^ 
that influences the exponential rate at which particles leave site x. In SEP a^ is 
exactly the rate of jumping from x, and in ZRP ct^ multiplies the rate r{ri{x)) that 
depends on the number r]{x) of particles currently occupying site x. The asymmetry 
pertains to the jump probabilities p{x,y), according to which a particle jumping 
from X chooses its new location y. We assume throughout that the kernel p{x, y) is 
translation invariant so that p{x, y) — p(0, y — x) = p{y — x). Asymmetric jumping 
means that there typically is a drift: 7 = J2^Pi^) ¥" 0- The assumption 7 7^ 
is not always necessary, but without it the limiting macroscopic conservation law 
becomes trivial. Some theorems require a stronger assumption of total asymmetry: 
the dimension d = 1, and all jumps proceed to the right: p{l) = 1. 

The disorder can also be attached to particles, so that individual particles carry 
their own randomly chosen jump rates. We do not explicitly consider such processes. 
One special case, the totally asymmetric simple exclusion process (TASEP) with 
particlewise disorder, is partially covered by our discussion. This is because the 
gaps between the exclusion particles with random rates can be regarded as the 
occupation numbers of a ZRP with random rates on the sites. This special case 
has been studied in the physics literature as a model for traffic. See Krug (1998), 
Krug and Ferrari (1996), and their references. 

The interesting phenomenon that appears in disordered particle systems is a 
phase transition where the process behaves differently at high and low densities. It 
occurs when the distribution of the random rate has a sufficiently thin tail at its 
left endpoint c > 0. Not much rigorous mathematical work exists on this phase 
transition. 

A brief overview of the paper: In Section 2 the disordered ZRP is described, 
together with two theorems. In Section 3 the same is done for the disordered 
SEP. Section 4 lists four open problems. Section 5 contains some proofs and some 
comments on proofs. In particular, we included in Section 5 a rigorous construction 
of a disordered ZRP in Z^, and a proof of the invariance of a certain family of 
product measures. The construction is based on the percolation approach of Harris 
(1972). 

The hope is that this paper would be at least partially accessible to the non- 
expert. This is the motivation for inclusion of the proofs in Section 5, which are 
often referred to but less often spelled out in the literature. For the same reason an 
attempt has been made to employ precise and complete notation. This may make 
the text somewhat heavy to follow at times, but the alternative is to risk confusing 
the reader who is not well-acquainted with disordered particle systems. Of course. 



such an outcome may be unavoidable in any case. 

Some familiarity with the subject of interacting particle systems is required for 
reading this paper. General references on particle systems are Durrett (1988, 1995), 
Griffeath (1979), and Liggett (1985). References on hydro dynamical limits are 
lectures by De Masi and Presutti (1991), the monograph of Spohn (1991), review 
papers by Ferrari (1994, 1996), and the soon-to-appear monograph of Kipnis and 
Landim. 

Here are some references that are closely related, but not directly on the topic of 
the paper: Hydrodynamic limits for asymmetric processes with inhomogeneous but 
not random rates have been proved by Landim (1996), Covert and Rezakhanlou 
(1997), and Bahadoran (1998). Koukkous (1996) and Giehs et al. (1998) have 
studied the symmetric ZRP with random rates. [In the symmetric case the jump 
probabilities satisfy p{x) = p{—x).] 

Notational remarks. Z_|_ = {0,1,2,3,...}. I^ and I{A} denote the indicator 
random variable of the event A. 5y is a delta function or a point mass at y, 
depending on the context: 5y{x) = l{x = y} for points x, and dx{A) = Ia{x) for 
sets A. 

2. The disordered asymmetric zero-range process 

First we describe a disordered ZRP on Z^ with bounded, monotone jump rates. 
Let {p{x) : X G Z*^} be a finite-range probability distribution, in other words 
p{x) > 0, '^p{x) = 1, and for some fixed finite set J\f C Z"^, p{x) = for a; ^ J\f. 
The rate of jumping from a site depends on the number of particles present through 
a function r : Z_|_ -^ [0, oo), about which we assume that 

(2.1) = r(0) < r(l) < r(2) < r(3) < r(4) < ■ ■ ■ 
and 

(2.2) r(oo) = lim r{k) < oo. 

k — >oo 

The disorder comes in the form of random deceleration factors a^ < 1 that depend 
on the sites x. Once a = {a^ '■ x E Z^) is picked, the dynamics operates as 
follows: If there are r]{x) > 1 particles at site x, then at exponential rate axr{r]{x)) 
a single particle jumps away from site x. The new location of this particle is y with 
probability p{y — x). This happens at each site x independently of what happens 
at other sites. 

For fixed rates ct, the generator of the process is 

(2.3) L^f{v)= E Piy-x)aMv{xm{v'''')-f{v)]- 



Here i] = {r]{x) : x G Z'^) is an element of the state space S = Z^ of the process, 
and 77^'^ is the configuration that results from the jump of a single particle from 
site X to site y: rj^'^ = rj + dy — dj;. Section 5 contains a construction of this process, 
based on a percolation argument of Harris (1972). Due to assumption (2.2) the 
process can be started from any configuration rj E S. We denote the process by 
r){t) — {r){Xj t) : X E Z'^), where t > is the time variable. 

The standing assumption is that a is an ergodic [c, l]-valued process for some 
constant c G (0, 1]. Let Q denote the distribution of the process a on the space 
A = [c, l]'^ . Fix c to be the left endpoint of the marginal distribution of ao, so 
that c is the largest number such that the process ex is [c, l]-valued. 

What makes the disordered ZRP tractable are invariant distributions that can 
be explicitly described. This description uses the same ideas as the process without 
disorder. If a^; = 1 (no disorder), among the extremal invariant distributions are the 
product measures fx^ on 5, indexed by a parameter i/j G [0,r(oo)), with marginals 



(2.4) ^^(^(^) = fc)=z(^)-i__ 



V' 



fc 



■r{k) 



X E Z'^ , k e Z+ . 



[See Andjel (1982).] For A; = the product in the denominator is interpreted as 1. 
Z{t(j) is the normalization factor, defined by 

For the disordered process, fix a choice a for the rates. For real numbers ip G 
[O,r(cxo)c), let u^ denote the product probability measure on S whose marginals 
vary from site to site, as given by 

(2.6) ^^{vi^) = k)=Z{^/a^)-'^,^^^y xGZ^fcGZ+. 

Proposition 1. For each choice of rates a E A and each value of ip E [O,r(oo)c), 
the probability distribution z/" is invariant for the process with generator L". 

The phase transition of the disordered ZRP is the following situation: If Q is such 
that very slow sites are sufficiently rare, then the family of invariant distributions 
{z/" : Lp E [O,r(oo)c)} does not cover the entire range of densities < p < 00. 
Instead, there is a critical density p* < 00 such that the equilibria z/" only exist for 
densities p E [0,p*]. To see this, set first 



Mi^i/j) is the density under fx^. It is a strictly increasing function from [0, r(oo)) 
onto [0, oo) and has an inverse function M~^ which we need to refer to below. For 
the disordered model the density p as a function of the parameter </? is defined by 
averaging over the random rates: 



(2.8) pM = E^ 



1 >^ kiip/ao) 



E;i^l-^«[M(./"o)]. 



Z{ip/ao){^r{l)---r{k) 



k=0 



Here E'^ denotes expectation over the distribution Q of ex, and the random variable 
inside the expectation is cto- For a fixed equilibrium z/" the density p{Lp) can be 
realized as a spatial average 

p((p) = lim -— ■ > ri(x) z/"-a.s., for Q-a.e. ex. 

^ AZ-zd A ^^-^ ^ 

By letting </? increase to its upper bound r(oo)c, (2.8) shows that the maximal 
density is 

(2.9) p* = E^ [M{r{oo)c/ao)]. 

This quantity may or may not be infinite, depending on the distribution Q. From 
(2.7) M(r(oo)) = lim^z-rj-Qo) M('i/') = oo, so in particular if Q{ao = c) > 0, then 
p* = oo. The interesting case with phase transition is the one where Q{dao) has a 
sufficiently thin tail as cto \ c, to make the integral in (2.9) finite. 

The function p : [O,r(oo)c) -^ [0, p*) is strictly increasing. Let / : [0,p*) -^ 
[O,r(oo)c) denote its inverse function. In other words, for p G [0,p*), /(p) is 
implicitly defined by 

(2.10) p = E^[M{f{p)/ao)]. 

Now we state a hydrodynamic limit for the disordered ZRP, due to Benjamini, 
Ferrari, and Landim (1996). For each choice of rates ck, there is a sequence of 
zero-range processes indexed by n, generated by L". P" denotes the probability 
measure on the probability space of the nth process {rjuix^t) : x G Z'^,t > 0}, 
n = 1,2,3,.... The theorem is a weak law of large numbers for the empirical 
measure defined by 

(2.11) TTnit) = n"^ Yl Vn{x,t)d^/r,, 

xeZ'i 
where S^ is a unit mass at the point x G R^. The assumptions are the following: 



(A.l) The transition probability p{x) satisfies this irreducibihty condition: for 
each Xj y E Z'^ there exists a finite sequence x = xq, . . . ^Xk = y such that p{xi+i — 
Xi) + p{xi — Xi+i) > for all i. 

(A. 2) There exists a bounded continuous function uq on R'^ such that ||wo||oo < 
M{r{oo)6) for some 9 < c, and for each a, the initial distribution of the process r]n 
is given by 

Pn{Vnix,0) =k) = HM~^uo(x/n)){'n{x) = k) . 

[Recall the definitions of /U^ and M{^) from (2.4) and (2.7).] 

(A. 3) The marginal distribution of cto is supported by a finite set: For some 



c 



ci < C2 < ■ ■ ■ < c^ < 1, (5(ao e {ci, . . . , Cm}) = 1. 



Assumption (A. 2) ensures that for some fixed </? G [0, r(oo)c) and all a, all the 
initial distributions are stochastically dominated by z/". This is true because, on the 
?7(a;)-marginal z/" = n^/a^, and this dominates Hr(oo)9 as long as v'/cix > r(oo)6', 
which in turn is true for all a^ € [c, 1] if </? > r{oo)9. 

Let 7 G R'^ be the mean drift under p{x): 7 = X]x6Zd^P(^)- Let w(x,t) on 
R'^ X [0, 00) be the unique entropy solution of the scalar conservation law 

(2.12) du/dt + -fV^[f{u)] = 0, u{x,0)=uo{x). 

Let Co(R'^) denote the space of compactly supported continuous functions on R''. 
This theorem was proved by Benjamini et al. (1996): 

Theorem 1. Under assumptions (A.1)-(A.3), the following holds for Q-a.e. a: 
For each t> 0, (p e Co{BJ^), and e > 0: 



nn{nt,(p)— I (p{x)u{x,t)dx 



> £ = 0. 



(2.13) lim P^ 

The integral against 7r^(t) is defined by 7r^(t, (/>) = n~'^'^^'rjn{x,t)(/){x/n). The 



/x 

>d 



statement is that, in the topology of Radon measures on R'^, iVn{nt) converges to 
u{x, t)dx in probability as n ^ 00. 

The shortcoming of this result is that it does not indicate what happens on the 
hydrodynamic scale if the process starts at density above critical, that is, uq{x) > p* 
for some or all x. In fact, assumption (A. 3) makes p* = 00, so there can be no 
phase transition under these hypotheses. 

Next we state a theorem that covers the hydrodynamics also in the high-density 
regime p > p* and admits more general initial distributions for the process. How- 
ever, we pay a serious price for this strengthening: The theorem is valid only for 



the most basic type of ZRP with rate function r{k) = l{k > 1}. Furthermore, we 
are restricted to totaUy asymmetric jumps in one dimension: d = 1 and p{l) = 1, 
so jumps happen only to the right on Z. FinaUy, we assume that the process of 
rates {ax : x G Z) is i.i.d. 

In this case the measures z/" are products of geometric distributions: 

(2.14) <(r?(x) = A;) = (1 - v?/«,)(v?/a.)'= , x e Z , k E Z+ . 

Now M{i(j) = t(j / {1 — t(j) , so the definition of the critical density becomes 



(tto - c) ^Q{dao). 

.1] 

From this formula it is plainly obvious how the tail of Q{dao) at ao = c+ determines 
whether p* < oo or not, that is, whether phase transition happens or not. 
For p G [0,p*) the flux function /(p) is deflned implicitly by the equation 



(2.15) 

If p* < oo, set 
(2.16) 



P = fip) («o-/(p)) ^Qidao). 

J[c,l] 



f{p) = c for p > p*. 



This makes / a nondecreasing and concave function on [0, oo). 

Again we assume we have a sequence of processes ?7n(t) and corresponding prob- 
ability measures P". The assumption on initial distributions is this: 

(A. 4) Suppose uq{x) is a nonnegative locally integrable function on R. Assume 
that this holds for Q-a.e. a.: For all (j) ^ C'o(R') and £ > 0, 



(2.17) 



hm P^ 

n — >oo 



TT, 



.(o,</>) 



(/){x)uo{x)dx 



R 



> e 



0. 



Previously this assumption was a consequence of assumption (A.l) so it was not 
stated explicitly. Let u{x, t) on R x [0, oo) be the unique entropy solution of 



(2.18) 



du/dt + df{u)/dx = 0, u{x,0) = uo{x) 



where / is the function deflned by (2.15)-(2.16). Then we have a theorem due to 
Seppalainen and Krug (1998): 



Theorem 2. Assume that d = 1, p{l) = 1, r{k) = l{k > 1}, and that Q is an 
i.i.d. distribution for a = (ax). Then under assumption (A. 4) the following holds 
for Q-a.e. ex: For each t > 0, (/> E Co(R); and e > 0: 



(2.19) 



lim P' 



nn{nt,(j)) 



^(x)w(x, t)dx 



R 



> e 



In the phase transition case this theorem has an interesting consequence: Sup- 
pose the initial profile satifies ^0(3:) > p* everywhere on R. Then since / is constant 
for this range of densities, it follows that u{x,t) = uo{x) for all t > 0. In other 
words, the profile does not change on the hydrodynamic scale. 

3. The disordered asymmetric exclusion process 

The disordered SEP is less well understood than the ZRP. The reason is that 
no invariant distributions have been found. Presently we can prove the existence 
of a hydrodynamic limit for the totally asymmetric SEP, in dimension one. But 
the fiux function / of the macroscopic equation (2.18) remains unknown. The 
theorem covers a more general totally asymmetric SEP, namely the so-called K- 
exclusion, where each site admits K particles instead of just one. The state space 
is 5 = {0, . . . jK}"^, and we write again r] = {r]{x) : x E Z) E S for the particle 
configurations. When the rates a have been chosen, the generator is 



(3.1) 



L'^fiv) = y , '^x'^{rj(x)>l,r](x+l)<K-l}[f{V^ 

xez 



x,x-\-l\ 



fiv)]- 



In other words, a jump occurs from site x to x -f- 1 at rate a^, provided site x is 
not empty and site x + 1 has less than K particles. Consider K fixed but arbitrary. 
As before, assume we have a sequence of totally asymmetric iiT-exclusion processes 
rin{t), with probability measures P" when the rates a. are fixed. As for Theorem 
2, we only assume that the initial distributions of the processes have a well-defined 
macroscopic profile: 

(A. 5) Suppose uq is a bounded measurable function on R such that < uo{x) < 
K. Assume (2.17) holds for all (j) G Co(R) and £ > 0, for Q-a.e. a. 



The theorem about the existence of the hydrodynamic limit is proved in Seppa- 
lainen (1998): 

Theorem 3. Fix a positive integer K and an i.i.d. distribution Q for ex. = (a^). 
Then there exists a concave function fx on [0, K] that depends on Q such that, 
if u{x,t) is the unique entropy solution of (2.18) with flux f = fx, then under 



assumption (A. 5), the limit in (2.19) holds for Q -a. e. a, for eacht > 0, ^ G Co(R); 
and £ > 0. 

The entropy solution to (2.18) with f = fx and uq bounded measurable can be 
constructed without explicit reference to the p.d.e: Pick a function Uq such that 
Uq = Uq a.e. on R. Let f^ be the concave conjugate of fK- 



U<p<-K 



Set 



(3.2) U{x,t) = snp{Uo{y) + tf*^{{x - y)/t)}. 

yen 

Finally, let u{x,t) = {d /dx)U{x,t), a derivative that is defined at least a.e. Lax 
(1957) discusses this construction of the entropy solution of a scalar conservation 
law. See also Section 3.4 in Evans (1998). The conjugate /^ has a probabilistic 
meaning in this context: It is the macroscopic shape of a growth model associated 
to the K-exclusion process [Seppalainen (1998)]. 

4. Open problems 

Problem 1. Extend Theorem 1 to describe hydrodynamic behavior also at high 
density p > p* . Or, equivalently, extend Theorem 2 to more general ZRP's. The 
proof of Theorem 2 in Seppalainen and Krug (1998) depends on special properties 
of the totally asymmetric ZRP with rate function r{k) = l{k > 1}. The proof of 
Benjamini et al. for Theorem 1 follows a strategy invented by Rezakhanlou (1991), 
and may be a better candidate for generalization. 

Problem 2. In the phase transition case, the hydrodynamic theorem reveals only 
trivial behavior at high density p > p* (the profile does not move, see Theorem 
2). It is expected that on a finer space scale one can see inhomogeneities develop, 
specifically, that the particles concentrate on exceptionally slow sites. No rigorous 
results exist to decribe these phenomena. Krug (1998) and Seppalainen and Krug 
(1998) discuss this in terms of the exclusion version of the totally asymmetric ZRP. 

Problem 3. An interesting open problem for the disordered ZRP concerns the 
weak convergence of the process: Fix the rates oc, and take the initial distribution of 
?7(0) spatially ergodic with density p. Does the distribution oir]{t) converge weakly, 
as t — » oo, to one of the equilibria z/"? If p > p* in the phase transition case, does 
r]{t) converge weakly to the equilibrium with density p*? For the standard ZRP 
results of this type were proved by Andjel et al. (1985). 

Problem 4- For the disordered TASEP, any new rigorous results would be wel- 
come. For example, is there a phase transition similar to the one for disordered 
ZRP? Does the fiux function fxip) have a fiat segment on an interval around 

p = K/21 
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5. Comments on the proofs 

5.1 Construction and invariant distributions. To construct the disordered 
ZRP with generator (2.3), start by giving each site x G Z*^ an independent rate 
r(oo) Poisson point process T^ on the time line (0, oo). [Here we make use of 
the assumption r(oo) < oo.] To the zth epoch of T^ attach two random objects: 
A random threshold Uf uniformly distributed on [0, r(oo)] and independent of 
everything else, and a destination site yf chosen with probability p(yf — x), again 
independently of everything else. Fix the rates o; = (a^ : x G Z*^) and an initial 
configuration r] = {r]{x) : x G Z'^). Informally speaking, from these ingredients the 
construction of the dynamics goes as follows: Suppose r is an epoch of T^ with 
threshold U and destination site y, and the dynamics rj{t) has been constructed for 
times < t < T. If U > axr{rj{T—)), do nothing. If U < axr{r]{T—)), move one 
particle from site x to site y; in other words, set 





rj{x,T) = rj{x,T-) -1, 




T]{y,T) = ?7(y,r-) + l. 


and 


r]{z,T) = r]{z,T-) for z^x.y 



Subsequently site x lies dormant until the next epoch of T^. However, the lattice 
Z'^ is infinite, so this induction never even begins because there is no first epoch 
among the point processes {T^ : x G Z*^}. 

To make the construction rigorous, we show that there exists a fixed time to > 
such that, starting with an arbitrary rj in the state space S, the evolution r]{t) can 
be computed for t G [0, to] • Since to is independent of ry, the construction can be 
repeated, starting with ?7(to), and extended to [0, 2to]. And so on, to arbitrarily 
large times. 

Given a fixed number to > and the Poisson processes {T^}, construct the 
following random graph, with vertex set Z^: Connect x and y with an edge if x — y 
or y — X lies in A/", and either T^ or T^ has an epoch in [0, to]. Recall that p{z) = 
for z outside Af. 

Lemma 5.1. For small enough to > 0, the random graph thus constructed has no 
infinite connected components, for almost every realization of {T^}. 

Before proving the lemma, let us observe how it solves the construction problem: 
All the sites y that could possibly influence the evolution at site x up to time to are 
connected to x in the random graph. Since x lies in a finite connected component 
C, the point process U^ec'^^ has only finitely many epochs in [0, to]. Consequently 
the evolution rj{Zjt) can be computed for 2 G C and t G [0,to] from the informal 
rule spelled out above, by considering the finitely many epochs in their temporal 
order. This procedure is repeated for all connected components. 
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Proof of Lemma 5.1. By translation invariance, it suffices to sfiow tfiat tfie origin 
is almost surely connected to only finitely many vertices. Let A/"* = A/" U (—A/") 
and K = | A/"* | , the number of vertices x such that either x or —x lies in A/". Call 
xo,xi, . . . ,Xn a self-avoiding path with n edges in the random graph if Xi 7^ Xj 
for i ^ j and there is an edge between Xi and Xi+i for each i. In particular, then 
Xi+i — Xi E Af* for each i. 

Let I ■ I be any norm on R*^, and R = max{|x| : x G Af*}. If the origin is 
connected to a vertex x with |a:| > L, there is a self-avoiding path with at least 
L/R edges starting at the origin. The probability that a self-avoiding path with 
2n — 1 edges starts at the origin is at most 

To see this, note first that the factor i^2n-i ^g ^^j^ upper bound on the number of such 
paths. If = xq, Xi, . . . , X2n-i is such a path, the n edges (xq, xi), (x2, 0:3), . . . , 
{x2n-2jX2n-i) a^e present independently of each other, and each with probability 
1 — e~2'^(oo)to i^g^^ igg^g^ Qj^g Qf q-x2i g^j^^ 7-x2«+i jjiust have an epoch in [0, to] , and each 
T^* has rate r(oo)]. Pick to small enough so that K'^(l — e~2'^(oo)to^ ^ i Then 
by Borel-Cantelli, self-avoiding paths from the origin have a finite upper bound on 
their length, almost surely. D 

This approach to the construction of a particle system is due to Harris (1972). 
Our presentation followed Durrett (1995). 

Let (O, JF, P) denote the probability space whose sample point to represents a 
realization of the Poisson processes {T^}, the random thresholds {U[}, and the 
destination sites {yf}. We constructed the random path ?](■) = {rj{x,t) : x G 
Z'^, t > 0} as a function of the initial state ry, the rates a, and a sample point to. 
Since the Poisson processes are Markovian, and the random choices of threshold 
U and destination state y are independent of everything, the process r]{-) is a 
time-homogeneous Markov process. Let V = ©([O, oo),iS) denote the space of 
right-continuous iS-valued functions on [0, 00) that have left limits at each point 
t G [0, 00). For fixed (ck, i]) the path r]{-) is a "D- valued random variable, and it has 
a probability distribution P"''? on V induced from the measure P on O. Write E"^'^ 
for the expectation under P"''?. 

Next we prove Proposition 1 about the invariant distributions. For this we 
restrict the process to a cube A^ = {— /c, . . . , /c}'^ C Z*^. The jump probabilities 
p{y — x) are then interpreted with periodic boundary conditions, and become 

(5.1) Pk{x, y) = 2^{p{z — x) : z = y + {2k + l)w for some w G Z*^} 
for x,y G A^. The finite- volume generator is 

(5.2) L'^f{v)= J2 Pk{x,y)aMv{xm{v^''')-f{v)]- 

x,yeAk 
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L" generates a Markov jump process with uniformly bounded rates on the countable 
state space Sk = 7i_^''. Existence of this process follows from standard textbook 
material. Write -B"'^ for expectations under distributions of the process restricted 
to Afc. Notice that if / is a cylinder function on S, L'^f = L'^f for all large enough 
k. 

Define dual transition probabilities by pl.{x,y) = pk{y,x) and p*{x) = p{—x). 
Then pi is obtained from p* exactly as pk from p according to (5.1). Let L"^'* and 
L"'* be the generators with kernels p^ and p* in place of pk and p. 

Lemma 5.2. Let ct E A and (f G [0, r(oo)c), and let u^ be the product probability 
measure defined by (2.14). Then for all bounded measurable functions f,g on Sk, 



-v 



■\gLtf)=v-{fLr9). 
Proof. Start by writing 

^^{9Lkf)= Yl Pk(x,y)a^iy^[r{r]{x))g{r])f{r]'' 

x,yeAk 

x,yeAk 
= A1-A2, 

where the last equality defines the quantities Ai and A2. 

Continue first with Ai: For any x,y let P denote the marginal distribution of 
fj = {r]{z) : z 7^ x,y). Abbreviate R{m) = r(l) ■ ■ -r^m). Then for a fixed x, the 
sum over y in Ai can be expressed as 

Y Pk{x, y)a^iy^ [r{r]{x))g{r])f{r]'''y)] 
yeAk 

= Y Pk{x,y)ax / J^idfj) Y Z{ip/aj:)~'^Z{ip/ay)~'^R{mx)~'^R{my)~'^ 

Tny>0 

{(p / aj:)'^'' {ip / ay)'^y r{mj:)g{fi,mj:,my) f {fi^ijij: - l,my + l) 
= Y Pk(^^y)^y / ^(^^) Y Z{ip/a^)~^Z{ip/ay)~^R{n^)~^R{ny)~^ 

yeAk n,j;>0 

ny >1 

{ip/a^)"'^{ip/ay)''yr{ny)g{fj, n^ + 1, % - l)f{f],n^, riy) 

= Y Pk{x,y)ayU^[r{r]{y))g{r]y''')f{r])]. 
yeAk 
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In the above calculation we first used the definition (2.14) of z/". The expectation is 
taken separately over the marginal distributions of fj, ri{x), and ri{y). In the second 
expression, rux and niy are summation variables that represent integration over 
the distributions of 'q{x) and r]{y). Because r{mx) = for rux = 0, we sum over 
"^x > 1. In the second equality we do a change of variable in the sum: n^ = nix — 1 
and Uy = ruy + 1. The last equality is just the definition (2.14) of z/" again. 
Now sum over x G Ayt and interchange the summation indices x,y to get 

^1= 5Z Pfc(2^,y)«x<[r(?7(x))(7(?7"'^)/(?7)]. 

x,y6Afe 

In A2 simply observe that YliyPk{x, y) = 1 = YlyPki^^ v) ^^^ write 
^2 = 5^ I XI Pk{x,y)\axi^^[r{r]{x))g{r])f{ri)\ 

x,yeAk 



Combining gives 



u-{gLtf)=A^-A2 = u^{fL'^'*g) 






and the lemma is proved. D 

Proof of Proposition 1. Let to > be the number chosen in Lemma 5.1. Since r]{t) 
is Markovian, invariance follows if we prove that, for any cylinder function / on 5, 



(5.3) [ E-'^[f{rj{t))]u^{dv)= j fdul 

Js Js 



for allt G [0,to]. 

For any rj E S, let tja,, denote its restriction to A^. Let ?7^(-) denote the process 
in Sk generated by L". (This of course is not the restriction of the process r]{-) to 
Afc.) Taking g = 1 in Lemma 5.2 shows that u^ is invariant for the process rj^{-), 
so we have 

(5.4) i^"'''' [^(^'(^))]<(^^) ^Is^'^''^ 

for all k large enough so that / can be regarded as a function on Sk- [Here again we 
rely on basic facts of continuous-time Markov chains on countable state spaces. The 
whole point is of course that the space S is not countable, so we need something 
more to pass from Lemma 5.2 to (5.3).] It remains to argue that 

(5.5) C'=[/(^'W)]-^"'n/(^W)] 
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boundedly, as /c ^ oo, for any fixed t G [0,to] and rj E S. For fixed (ck,?]), we can 
construct the processes ?]{■) and r)'^{-) on the same probabihty space (0,JF, P). We 
only need to interpret the destination sites yf "modulo A^" for the process r)^{-). 
Then both f{rf{t)) and f{ri{t)) are functions on Q, and both expectations in (5.5) 
are integrals over (0,JF, P). Since / is bounded, it suffices to show that 

(5.6) f{Ti\t)) ^ /(r7(t)) as fc ^ oo, 

almost surely on (n,jF, P). 

Pick /cq large enough so that firf) is completely determined by {r]{x) : x G 
Afcg ) . Let Cjto be the random set of vertices connected to A^^ in the random graph 
discussed in Lemma 5.f. By that lemma, Ck^ is almost surely finite. But then 
f{vi^{t)) = f{T]{t)) as soon as k is large enough so that C^y C A^, because then all 
the transitions that affect the value of / are performed identically for r][-) and r]^{-) 
throughout the time interval [0, to]- This proves (5.6) and thereby (5.5), and then 
(5.4) turns into (5.3) as /c ^ oo. D 

5.2 The hydrodynamic limits. Theorem 1 is Theorem 3.1 from Benjamini et 
al. (1996). It is proved by deriving a microscopic version of Kruzkov's entropy 
inequalities that characterize the entropy solution of a conservation law. This idea 
for proving hydrodynamic limits of asymmetric processes is due to Rezakhanlou 
(1991). 

Theorem 2 follows from Theorem 2 in Seppalainen and Krug (1998). The discus- 
sion in this paper is formulated for a totally asymmetric exclusion process where the 
random rates are attached to the particles. To obtain results for the ZRP with site 
disorder, interpret the occupation numbers (ry^) as the gaps (= number of empty 
sites) between successive exclusion particles. The proof uses a special construction 
that works for the totally asymmetric ZRP with rate function r{k) = l{k > 1}. 
Whether the technique can be somehow generalized to deduce results for other 
choices of r{k) is presently unclear. 

Theorem 3 is proved in Seppalainen (1998). The approach is similar to the one 
in Seppalainen and Krug (1998). It involves coupling the exclusion process with 
a countable collection of realizations of the same process but with a simple initial 
configuration. This coupling has the property that a microscopic version of the 
Lax-Oleinik formula (3.2) holds almost surely. 

Acknowledgment. 1 thank IME at the University of Sao Paulo for fruitful hos- 
pitality, and Pablo Ferrari, Joachim Krug, Claudio Landim, Gunter Schiitz, and 
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